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Abstract 

A scenario of the big-bang nucleosynthesis is analyzed within the minimal supersymmetric stan- 
dard model, which is consistent with a stau-neutralino coannihilation scenario to explain the relic 
abundance of dark matter. We find that we can account for the possible discrepancy of the abun- 
dance of Li between the observation and the prediction of the big-bang nucleosynthesis by taking 
the mass of the neutralino as 300 GeV and the mass difference between the stau and the neutralino 
as (100- 120) MeV. We can therefore simultaneously explain the abundance of the dark matter 
and that of 7 Li by these values of parameters. The lifetime of staus in this scenario is predicted to 
be 0(100-1000) sec. 
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The supersymmetric models are attractive candidates of the theory beyond the standard 
model. While no experiments so far have found any evidence of the supersymmetry , the 
Large Hadron Collider is expected to find its first signal in the near future. On the other 
hand, the analysis of the cosmological implications of the supersymmetry is an approach 
complementary to the direct search. The lightest supersymmetric particle (LSP) is stabilized 
by the R parity and naturally qualifies as the cosmological dark matter. A possible candidate 
of the LSP, and hence of the dark matter, is the neutralino x°- The neutralinos as the LSP 
with their mass of 0(100 GeV) can be responsible for the present abundance of the dark 
matter when the mass of the next-lightest supersymmetric particle (NLSP) is close to that 
of the LSP's and allow them to coannihilate with each other in the early universe. 

We put this coannihilation scenario into the perspective of the big-bang nucleosynthesis 
(BBN). The recent results from the Wilkinson Microwave Anisotropy Probe experiment 
combined with the standard BBN scenario, suggest twice or thrice as much abundance of 7 Li 
as suited fr0 ,„ the option of m e t al-p„c, halo stars QflQ. This discrepancy may 
imply that the 7 Li nuclei were destructed in the BBN era through processes of physics beyond 
the standard model although there might still be a possible astrophysical process to deplete 
7 Li uniformly . We introduced in Ref . [f| a scenario in which an exotic negatively- charged 
massive particle form a bound state with a nucleus and therethrough initiate the destruction. 
There we analyzed the Minimal Supersymmetric Standard Model with the coannihilation 
scenario where the NLSP is the stau f. Staus serve as charged massive particles that trigger 
the destruction of 7 Li through the interaction 

AC 

Ant = t*x°(<7lPl + <?rPr)t + — |i/ T fP LT jj ad 
+ —^(l^ PjjUl )(u TjlM P LT ) + H.c, 

where Gp = 1.166 x 10~ 5 GeV~ 2 is the Fermi constant, Pl and Pr are the chiral projec- 
tion operators, I G {e,/i}, g L and g R are the coupling constants, and J^ ad is the hadron 
current. The stau f in Eq. ([1]) is the mass eigenstate, which is given by the linear combi- 
nation of the superpartner of the left-handed tau tl and that of the right-handed tau tr as 
f = tl cos r + fRe~ 17r sin 6 T , where 6 T and 7 T are the left-right mixing angle and CP- violating 
phase, respectively. The formation of a stau- 7 Be bound state f + 7 Be — > (f 7 Be) + 7 is imme- 
diately followed by an internal conversion process (f 7 Be) — * x° + v t + 7 Li, and subsequent 
spallation of 7 Li by the energetic protons in the background. We assumed in Ref. ||3] the 



rapid formation of the stau-nucleus bound state and ignored the effect of the expansion of 
the Universe, as the use of the Saha equation implies. 

In the present paper, we improve the previous analysis by considering the expansion 
effect of the Universe. The Boltzmann equation is employed instead of the Saha equation to 
estimate the stau-nucleus bound states. We also include the resonant formation of the bound 
state pointed out in Ref. Q]. We thereby show that the LSP and NLSP, both with mass 
of 0(100 GeV), can account for the problem of the dark matter and that of the abundance 
of 7 Li. The relevant parameters in considering our BBN scenario are the mass difference 
between staus and neutralinos 8m = — m^o, where rrif and m^o are the masses of 
staus and neutralinos, respectively, and the yield value of the staus at the freeze-out time 
Y^fo = riffs | Freeze Out, where rif and s are the densities of the number of staus and the 
entropy, respectively. Other parameters are fixed throughout this paper to 9 T = 7r/3, 7 r = 0, 
and m^o = 300 GeV. By varying the values of 5m and Yf^o, we search for the parameter 
region that can account for the present abundance of 7 Li. Staus play a major role in the 
BBN when Sm is small so that staus become longevous enough to survive until the BBN 



era. The lifetime of staus indeed becomes 100 sec or longer when 5m < 100 MeV 

We show in Fig. [1] the evolutions of the bound ratios of He, 7 Li, and 7 Be, where we 

define the bound ratio by the number density of a nucleus that forms a bound state with 

a stau, divided by the total number density of that nucleus. We trace the evolution of the 

number density of the stau-nucleus bound states by the Boltzmann equation, using the cross 

sections shown in Ref. 9|. 

The yield value of staus at the time of the formation of the bound state with nuclei £bf, 

which we denote by Y f ^p, is changed from 10~ 10 to 10~ 16 in each figure. It is related with 

Yf using the lifetime of stau Tf as 

n,BF = n,Foe-' BF/T? . (2) 



The bound ratio of 4 He shown in Fig 
3 Li due to the catalyzed fusion process 
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a) is crucial to estimate the creation rate of 
ll| (f 4 He) + D -> 6 Li + f, while that of 7 Be 



shown in Fig. QJb) is necessary to evaluate the reduction rate of the 7 Li. Since the present 
7 Li originates from the primordial 7 Be, the abundance of 7 Li is reduced as follows. We first 
convert 7 Be into 7 Li by an internal conversion process (f 7 Be) — > x°+i/ r + 7 Li, and successively 
destruct the daughter 7 Li by either a collision with a background proton or a subsequent 
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FIG. 1: Evolutions of the bound ratio of the nuclei 4 He, 7 Be, and Li. We vary the abundance of 
the stau at the time of the formation of the bound state from 10 _ 10 to 10 _ 16 in each figure. In 
Fig. Ufa), we also plotted by dotted lines corresponding curves predicted using the Saha equation 
for reference. 

internal conversion (f 7 Li) — > x° + v T + 7 He. The 7 Be is efficiently reduced if its bound ratio 
plotted in Fig. QJb) is of 0(1). The successive destruction of 7 Li by internal conversion is 
effective when its bound ratio plotted in Fig. [Tfc) is also of 0(1). We find in Figs. Hfb) and 
(He) that both bound ratios of 7 Be and 7 Li are of 0(1) when Yf^ FO > (10~ 13 - 10 -12 ). 

The parameter region that can solve the 7 Li problem is numerically calculated in the 
(8m, Yf t po) plane and presented in Figure [21 in which Fig. [2( a) does not include effects 
of the resonant formation and photo-dissociation processes of the bound state pointed out 
in Ref. [3], while Fig. 12(b) includes these effects. The white region is the parameter space, 
which is consistent with all the observational abundance including that of 7 Li/H. The region 
enclosed by dashed lines is excluded by the observational abundance of 6 Li/ 7 Li [jjj], and the 
one enclosed by solid lines are allowed by those of 7 Li/H 3|. The thick dotted line is given 
by the upper bound of the yield value of dark matter 

0.110 J Vl0 2 GeV, 



taking f2 DM /i 2 = 0.1099 + 0.0124 (upper bound of 95% confidence level) [l| and m DM 
This line gives the upper bound of I^fcm since the supersymmetric particles after their 
freeze-out consist of not only staus but neutralinos as well in our scenario. 

The allowed region shown in Fig. [2] lies at 8m ~ (100-120) MeV, which is tiny compared 
with rriyA) = 300 GeV. These values of parameters allow the coannihilation between neu- 
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FIG. 2: Allowed region in bm-Yf plane. The resonant capture and the photo dissociation are 
(a) not included, and (b) included. The white region is the parameter space, which is consistent 
with all the observational abundance including that of 7 Li/H. The region enclosed by dashed lines 

n 

is excluded by the observational abundance of 6 Li/'Li [12|, and the one enclosed by solid lines are 
allowed by those of 7 Li/H j|. The thick dotted line is given by the upper bound of the yield value 
of dark matter. This line gives the upper bound of If fo- 



tralinos and staus, and thus can account also for the abundance of the dark matter. We 
therefore find that the values of m^o = 300 GeV and 5m ~ 100 MeV can simultaneously 
explain the abundance of dark matter and of 7 Li. 

We compare the Figs. |2^a) and [2(b) to find that the allowed region is shifted downward in 
Fig. [2](b). Of the two processes included in Fig. [2(b), the resonant formation process makes 
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the bound ratio larger, the value of Yf smaller, and push the allowed region downward in 
the figure. On the other hand, the photo dissociation process makes the bound ratio smaller 
through the destruction of the bound state, makes the value of Yf larger, and push the 
allowed region upward. We thus find that the resonant formation of the bound state is 
relevant while the photo dissociation is inconsequential. 

The qualitative feature of the allowed region is explained from the following physical 
consideration. First, we note that V^fo ^ (10 -13 - 10~ 12 ) is required so that a sufficient 
number of bound state (f 7 Be) is formed to destruct 7 Be by the internal conversion into 
7 Li. The daughter 7 Li is broken either by an energetic proton or by the internal conversion 
(f 7 Li) — > x° + v T + 7 He, and consequently 7 Li/H is reduced. Bearing this physical situation 
in mind, we consider parameter regions in detail. 

1. 5m > 120 MeV. 

Since the staus decay before they form a bound state with 7 Be, the value of If bf is 
much lower than 10~ 13 and hence the abundance of neither 7 Be nor 7 Li is reduced. 
Therefore this parameter region is excluded. 

2. 100 MeV < 5m < 120 MeV. 

The staus are just decaying at the formation time of the bound state. The necessary 
condition of V^bf ~ 10 -13 can still be retained even in a case where the value of I^fo 
is sufficiently large. The allowed region in this area of 5m thus bends upward. In this 
region, a daughter 7 Li from the internal conversion of (f 7 Be) is broken mainly by an 
energetic proton. 

3. n, FO < 10- 13 . 

In this case V^ bf is necessarily less than 10 -13 , and the bound ratio of 7 Li and 7 Be 
are much less than 0(1) as seen in Fig. [TJ Therefore, the final abundance of 7 Li is not 
reduced sufficiently. This parameter region is thus excluded. 

4. Y f)FO > 10~ 12 and 5m < 100 MeV. 

In this region Yf.BF — ^t,fo > 10~ 12 and hence the bound ratio of 7 Be is 1 (see Fig. [I]). 
It means that 7 Be and consequently 7 Li are destructed too much. Hence, the upper-left 
region is excluded. 

5. 5m < 100 MeV and Y f>F0 > 10~ 15 . 
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In this region, the stau acquires the long lifetime enough to form a bound state (f 4 He) . 
Then the catalyzed fusion process (f 4 He) + D — ► 6 Li + f leads to the overproduction 
of 6 Li and to the disagreement to the observational limit. Therefore, this parameter 



region is excluded, which is consistent with calculations by Ref. 
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Excluding all the parameter regions described above, we obtain a small allowed region of 
m^o ~ rrif ~ 300 GeV and 5m = (100-120) MeV as presented in Fig. [2J and these values 
are at the same time consistent to the coannihilation scenario of the dark matter. 

We obtained a strict constraint on the mass of the neutralinos and staus by improving 
an analysis of a solution to the overproduction problem of 7 Li or 7 Be through the internal 
conversion in stau- nucleus bound states, (f 7 Be) — > x° + v t + 7 Li and (f 7 Li) — > x° + v r + 7 He, 
given in Ref. {(J. We included the resonant capture process of 7 Be and photo-dissociation 
process pointed out in Ref. 3]. We also took into account the expansion of the Universe by 
an explicit use of the Boltzmann equation instead of the Saha equation to obtain a more 
accurate number of the stau-nucleus bound states. By varying the yield value of the stau 
at its freeze-out time, we found that most of 7 Li and 7 Be nuclei form a bound state with a 
stau for Yf t BF ^ (10~ 12 - lO" 1 ^). Taking the values of m^o = 300 GeV, 9 T = ir/3, 7 T = 0, and 
7] = (6.225 ± 0.170) x 10~ 10 [lj, we compared the primordial abundances with and without 
the resonant capture and/or photo dissociation, and found that the resonant capture process 
is relevant while the photo dissociation process of the bound state is inconsequential. We 
obtained a parameter region consistent with the observed abundance of 7 Li within Yf -po — 
(7-10) x 10~ 13 and 5m = (100- 120) MeV. The region of 5m < 100 MeV is excluded 
due to the overproduction of 6 Li by the catalyzed fusion. Furthermore, the parameter 
region obtained in this paper lies in the coannihilation region, which can explain the relic 
abundance of dark matter. Therefore, the stau with mf ~ 300 GeV and 5m ~ 100 MeV 
can simultaneously solve the problems on the relic abundance of the light elements and the 
dark matter. As shown in Ref. jg], the stau with m f = 300 GeV and 5m = 100 MeV has the 
lifetime of O(100- 1000) sec. It is very possible that Large Hadron Collider will find some 
staus with a very long lifetime fl^ |. 

We need further improvement on our analysis to obtain a more precise result of the 
mass and the mass difference. We have to derive Y? as a function of the parameters in the 
Lagrangian, although we regarded Yf as a free input parameter in this paper. Then, we can 
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determine the allowed region of 8m and m^o more precisely by varying other parameters 
such as 9 T and 7 T . We leave this for our future work. 
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